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对 不 可 压缩 两 相 闯 流 在 突变 界面 处 连接 条 件 的 辨析 与 重建 
B/N TE) 朱 苏 阳 ! EEX! 梁 保 升 ? 邓 鹏 13 机 春生 ! 
(1. 西 南 石油 大 学 油气 地 质 及 开发 工程 国家 重点 实验 室 ; 610500; 2. University of Texas at Austin;3. 加 拿 
大 卡尔 加 里 大 学 ) 
【摘要 】 真 实 油气 藏 儿 乎 都 有 突变 界面 ， 建 立 相 应 的 渗流 数学 模型 ， 其 关键 是 突变 界面 处 的 连接 条 件 。 已 
有 的 渗流 理论 已 确定 根据 质量 守恒 定律 ， 在 突变 界面 处 各 相 渗流 速度 和 压力 都 连续 《CPVCM) 。 但 通过 
调研 文献 发 现 : ORG CPVCM 推断 的 突变 界面 两 侧 的 流体 饱和 度 分 布 存 在 实际 的 反例 ; OERE 
CPVCM 与 Rankine-Hugoniot 跃 变 条 件 矛 盾 ，@@) 由 它 推导 出 的 数值 计算 格式 被 油 藏 数值 模拟 证 明 会 出 现 非 物 
理 结果 而 被 弃 用 。 继 而 回溯 了 文献 根据 质量 守恒 方程 推导 CPVCM 的 过 程 ， 发 现 这 些 推导 过 程 最 终 都 出 现 
指 代 错 位 ， 即 用 界面 上 同一 质点 的 流速 代表 了 界面 两 侧 不 同 质点 的 流速 ， 实 质 上 预 置 了 渗流 速度 连续 ， 属 
于 自 证 。 接 着 以 不 可 压缩 流体 为 例 ， 同 样 根据 质量 守恒 定律 和 两 相 流 达 西 公式 推导 出 各 相 流 速 和 压力 都 不 
连续 的 界面 条 件 (JPVCM) ， 从 而 建议 不 可 压缩 两 相 渗流 应 将 上 游 权 函 数 和 总 流速 相等 作为 突变 界面 处 的 
连接 条 件 。 
关键 词 : 突变 界面 ， 界 面条 件 ; ”两 相 流 ， 多 和 孔 介质 ; 质量 守恒 
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Analyzing and Rebuilding Sharp Interface Conditions for Incompressible Two-Phase Flow in Porous Media 
Peng,Xiaolong!;Zhu Shuyang!; Wang Chanwen!;Liang Baosheng’; Deng Peng!; Jia Chun-Sheng! 
(1. Southwest Petroleum University. State Key of Laboratory of Oil and Gas Reservoir Geology and exploration; 2. 
University of Texas at Austin, U.S. ; 3. University of Calgary, Canada) 


Almost all the oil or gas reservoirs have sharp interfaces. The key to establish the corresponding mathematical 
flow model is to find the reasonable phase pressures and velocities connection conditions at the sharp interfaces. In 
the existing seepage mechanics theory, it is generally recognized that the velocity and pressure of each phase of fluid 
are continuous at the sharp interface (named CPVCM). However, we can find that CPVCM contradicts some other 
multiphase seepage phenomena or theories: (1) The saturation distribution of fluid in the real world does not always 
obey the fluids distribution rules that CPVCM required; (2) Near the flood front, CPVCM conflicts with Rankine- 
Hugoniot Interface Conditions; (3) It has been proved that the formula of phase fluxes across the sharp interface 
which are derived from CPVCM, harmonic average of transmissibility, possible give non-physical results in some 
cases, so that they have been replace by single points upstream weighting method(SPU) and other high order method, 
such as TVD,ENO,WENO et al. . Therefore, we retraced the way that how to derive the CPVCM based on the mass 
conservation law in the early lectures, and find that the two particles on both sides of the interface which mentioned 
in the interface condition has been misplaced on the same interface, therefore, the proof process is essentially 
equivalent to preset the continuity of velocity and pressures of each phase and prove they are continuous. Such proof 
is invalid self-proof. Then, taking the incompressible two-phase fluids flow in the porous media as an example, the 
interface condition with discontinuous velocity and pressure of each phase (JPVCM) is obtained according to the 
same laws of mass conservation and two-phase Darcy’s law. Finally, it is suggested that, as for the incompressible 
two-phase flow in porous media, the jump conditions should be (1) total velocity of the two phases are equal, (2) the 


velocity of each phase at the sharp interface adopt the upstream ones. 


Keywords: Sharp Interface; Interface Conditions; Multiphase; Porous Media; Mass Conservation 
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油气 藏 突变 界面 是 指 岩 石 属性 、 流 体 分 布 或 性 质 在 某 些 空间 曲面 邻近 的 狭窄 区 域 急剧 变化 ， 该 空间 曲 
就 是 突变 界面 。 突 变 界 在 在 真实 油气 藏 中 广泛 存在 : 储 层 界面 、 沉 积 相 界 面 、 断 层 、 裂 颖 - 基 岩 等 跨 介 质 
类 型 跨 介 质 尺 度 界面 、 压 裂 颖 /井筒 - 储 层 界 面 ;油气 藏 原始 气 - 油 -水 界面 、 注 水 / 气 / 聚 合 物 / 示 踪 剂 的 驱 蔡 前 
缘 、 边 底 水 或 气 顶 锥 进 界面 、 流 体 相 变 界 面 。 总 之 ， 真 实 油气 藏 几乎 都 有 突变 界面 ， 指 导 油 气 藏 科学 开发 ， 
油 藏 工程 理论 必须 基于 反映 突变 界面 特征 的 渗流 数学 模型 ， 而 建立 合理 的 突变 界面 条 件 是 其 中 的 关键 。 
在 渗流 力学 领域 自 多 相 渗流 力学 建立 不 久 凹 ， 突 变 界面 条 件 就 已 建立 PC71， 即 : 各 相 流体 的 渗流 速度 和 
流体 压力 在 突变 界面 处 连续 〈 CPVCM) 已 *9。 数 在 学 上 ， 有 突变 界面 的 渗流 问题 属于 典型 的 Neumann 问 
题 00， 按 理 只 要 建立 CPVCM 结合 连续 介质 场 的 渗流 理论 就 足以 描述 考虑 了 突变 界面 的 渗流 问题 。 但 是 ， 
技术 调研 和 文献 调研 发 现 CPVCM 基本 没有 有 具体 应 用 于 解决 突变 界面 相关 的 油 藏 工程 问题 ， 即 使 这 类 问题 
得 到 解决 也 与 CPVCM 没有 密切 关系 ， 例 如 : 构建 油 藏 数值 模拟 差分 格式 、 预 测 驱 蔡 前 缘 位 置 、 计 算 裂 颖 - 
EARME MARIES, FH CPVCM 还 与 诸多 渗流 现象 或 理论 相 矛 盾 。 这 不 禁令 人 担忧 CPVCM 
的 合理 性 ， 反 思 它 的 科学 性 。 为 此 我 们 对 其 开展 了 调研 工作 ， 论 文 内 容 安排 如 下 : 首先 在 第 2 节 从 多 个 案 
例 剖析 CPVCM 存在 的 问题 ;然后 在 第 3 节 回 蛮 传 统 渗流 力学 为 什么 由 质量 守恒 定律 推导 出 了 CPVCM; 第 
—_ 4 节 以 不 可 压缩 两 相 流 为 例 ， 仍 按 质量 守恒 规律 和 达 西 公式 ， 但 得 到 了 不 同 于 CPVCM 的 突变 界面 条 件 〈 记 
a 为 JPVCM); 最 后 探讨 了 多 相 渗 流 方程 在 不 连续 界面 处 本 应 该 的 连接 条 件 。 

本 ”2. 现 有 突变 界面 条 件 存 在 的 问题 

= 突变 界面 作为 油 茂 的 内 部 界面 将 油 藏 分 割 为 多 个 子 区 域 。 描 述 油 藏 渗流 规律 的 数学 模型 包括 两 部 分 : 
(1) 描述 子 区 域 渗流 规律 的 数学 微分 方程 ， 因 子 区 域 连续 ， 所 以 它们 就 是 常规 连续 场 的 渗流 数学 微分 方 
程 ，(2) 突变 界面 的 连接 条 件 是 将 相 邻 介质 的 热力 学 性 质 联系 起 来 的 界面 条 件 ， 这 种 关系 通常 又 被 称 为 
“跳跃 条 件 ” 或 “突变 条 件 ”(Jump Conditions) |), 


图 1 突变 界面 将 油 藏 分 割 为 两 个 连续 的 子 区 域 及 渗流 数学 模型 
Fig. 1 the sharp interface divides the reservoir into two continuous sub regions 
( 工 表示 突变 界面 ，T TA 无限 靠近 昌平 行 荆 的 界面 E、F 分 别处 于 突变 界面 两 侧 的 两 点 ;2 Q, 
指 被 突变 界面 分 隔 开 的 两 个 内 部 连续 的 子 区 域 ) 
根据 文献 两 种 突变 界面 处 有 界面 连续 条 件 B-7; 


[o (v -w)|N=0 ef{n,w} (1) 


式 中 pp 表示 由 根据 微观 尺度 下 流体 分 布 计算 得 到 的 平均 密度 ; Vw 表示 上 相 流 体 的 流速 ，w 表示 间断 界面 的 
移动 速度 ;N 表示 界面 的 单位 法 向 矢量 。 下 标 n 表示 非 润 湿 相 ; w 表示 润 湿 相 ， 对 于 不 可 压缩 流体 ， 密 度 
为 常数 ， 公 式 O) 所 示 界 面条 件 可 化 简 : 

[(v. -w) |N =0 le{n,w} (2) 


对 于 空间 固定 界面 W=0 所 以 有 : [v |N =0 or v =v; le {n, w} (3) 
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即 : 各 相 流体 的 渗流 速度 和 流体 压力 在 突变 界面 处 连续 ， 为 方便 表述 将 公式 〈1) - (4) 记 为 CPVCM。 由 


CPVCM 和 达 西 公式 ， 考 虑 毛管 力 和 重力 作用 ， 可 得 到 压力 和 毛管 力 的 连接 条 件 225; 


vi =v Pi =P P: = pb Let{n,w} (4) 


CPVCM *U Axe ta tt PE A Ae, TEA Wa Ti eS OS, JL 
ZERI SHY ACE SEP 225], (AY ACE CPVCM 与 一 些 渗流 现象 和 其 它 的 一 些 渗流 理论 存在 不 可 调和 的 


矛盾 ， 导 致 我 们 对 其 质疑 。 


为 了 突出 问题 的 核心 ， 论 文 以 不 可 压缩 两 相 渗流 为 例 进行 说 明 ， 渗 流 运动 方程 服从 达 西 定律 ， 考 虑 习 


力 和 毛管 力 〈 见 公式 〈5)): 
v, =-KA,(Vp,-p,.gVD) Ce{n,w} (5) 


pan 


HP D=(0,0,z) 。 根 据 n，w 两 相 流体 的 质量 守恒 方程 ，Wp,9,+V.w, =0,00,S,+V-v, =0 得 到 


V-(v,+v,)=0 (6) 
并 定义 : 


v, =V; +V,» A =k, l Ul, A =A, PAS =A, lf, =1-f, 


AP = Py PrP = fo Po + fapa =24,A, 14, 4D = D-Day 

根据 等 式 (6) 可 知 Y 在 渗流 空间 为 无 散场 ， 为 连续 场 ， 则 达 西 公式 也 可 用 一 了 ,公式 改写 为 : 
Vide Gh. Va ERG, (8) 
其 ! GP =0.5K4,,, (Vp, +ApgVD) (9) 


2.1 质疑 1: 在 均 质 储 层 中 CPVCM 条 件 与 流体 饱和 度 场 的 矛盾 


不 考虑 毛管 力 的 一 维 两 相 不 可 压缩 流体 渗流 问题 即 为 Buckley-Leverete (B-L) 模型 [tl， 该 模型 能 够 得 


到 的 饱和 度 和 流速 场 的 解析 解 ， 解 析 解 显示 在 驱 蔡 前 缘 流 体 处 含水 饱和 度 和 水 相 流 速 不 连续 ， 如 图 


显然 不 符合 CPVCM 对 各 相 流 速 在 任何 界面 处 连续 的 要 求 。 为 了 解释 该 矛盾 ， 一 些 学 者 提出 05 27); 


2 所 示 ， 
真实 的 


多 和 孔 介质 总 有 毛管 力 ， 如 果 考 虑 毛管 力 ， 驱 替 前 缘 出 现 两 种 流体 的 过 渡 带 呈 连 续 分 布 ，B-L 模型 只 是 
对 于 是 


化 模型 。 我 们 认为 通过 和 否定 B-L 模型 来 支持 CPVCM 的 理由 非常 牵强 。 首 先 CPVCM 的 建 模 过 程 ， 
否 考虑 毛管 力 没 有 限定 ， 得 到 CPVCM 的 结论 与 是 否 考虑 毛管 力 无 关 ， 重 复 CPVCM 的 推理 过 程 ， 


理想 


即使 不 


考虑 毛管 力 同样 得 出 CPVCM， 所 以 B-L 模型 应 当 作为 CPVCM 的 反例 ; 其 次 ，B-L 模型 和 CPVCM 基于 


相同 的 物理 模型 ， 基 于 相同 的 质量 守恒 方程 、 运 动 方程 、 状 态 方程 以 及 定 解 条 件 ， 理 应 得 到 相同 的 数学 解 ， 
但 两 者 却 得 出 截然 不 同 的 突变 界面 条 件 ， 因 为 B-L 通 过 引入 也 三 多 十 y, ， 将 不 连续 场 问题 转换 为 连续 场 问 


题 ， 所 以 B 工 更 有 可 信和 度 。 我 们 提供 两 个 即使 有 毛管 力 ， 在 也 存在 流体 饱和 度 间 断 的 直观 的 反例 : 


反例 1. 如 图 3 所 示 ， 两 块 性 质 相 同 岩 石 两 种 流体 按 不 同 饱 和 度 充满 孔 除 ， 且 平衡 分 布 ， 然 后 足够 快速 
地 将 两 块 岩 石上 下 又 置 ， 在 重力 作用 和 毛管 力作 用 下 流体 将 重新 分 布 建 立新 的 平衡 ， 但 新 平衡 的 建立 过 程 


需要 时 间 ， 所 以 在 两 块 岩石 接触 的 瞬时 ， 界 面 处 流体 饱和 度 的 分 布 是 不 连续 的 。 
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Xf 


图 2 Buckley -Leverate 模型 水 驱 油 前 缘 饱 
和 度 和 水 相 流 速 示 意图 
Fig.2 Schematic diagram of water flooding 


front saturation and water phase velocity of 
Buckley leverate model 


反例 2. 如 图 4 所 示 ， 我 们 设 


different water saturation 
计 了 一 个 利用 均匀 球体 结合 
在 于 均 球 最 小 喉 道 半径 相对 孔 腹 3 
和 上 端 ， 分 别 取 毛管 力 | 


两 岩 块 瞬间 上 


下 拼接 


图 3 两 块 合 水 饱和 度 不 同 的 相同 宕 块 瞬间 
下 接触 示意 图 


Fig. 3 Schematic diagram of instantaneous up- 


T 


down contact of two same type of rocks with 


ae AA 


Fl PE Wit AS AT EA SR FTA SS 
# 径 较 大 ， 导 致 毛管 力 曲 线 出 现 较 长 水 平 段 。 在 重力 作 月 
F 段 的 最 大 值 和 最 小 值 ，》 


1 线 水 了 


日 下 油水 界面 下 端 
水 界面 就 成 为 饱和 度 的 突变 界面 。 


cowater =a oil @ solid particle 
利用 均匀 球体 设计 多 孔 介 质 


0s” 


: sr x, 
油水 分 布 示意 图 油水 分 布 示意 图 
图 4 一 种 制作 流体 饱和 度 呈 间断 分 布 的 实验 方法 
一 步 : 在 t=0 时 刻 让 整个 ; 


Fig. 4 an experimental method for making discontinuous fluid saturation field. 


水 系统 达到 垂 向 平衡 状态 ， 根 据 公 式 8) 和 CO) A: 


-和 (SBpe(S2)+Apg |=0 
= (10) 
-74 (s#)[ a.p. (s7 )+Apg]=0 
第 二 步 : 在 垂 向 方向 给 多 孔 介质 区 域 施 加 压 差 ， 让 流体 自 下 往 上 流动 ， 总 的 流速 w >0 ， 在 开始 流动 的 瞬 
间 ， 流 体 饱和 度 与 六 0 时 刻 相 同 ， 毛 管 力 与 重力 作用 仍然 平衡 ， 所 以 根据 公式 (8) 和 (9) 有 : 
v, =v = fo (S? )v, -0.5K Zn (Se)| 6.p. (S2")+ Avg |= F,(52"), (11) 
vi =v? = fy (S2 )v, -0.5K 3n (9%) 8,7, (S7 )+ Ape |= fo (52 )y, (12) 
由 于 8w eS. AS, (S,,) BE Sw AiR, BELLA: 
vv: 同 理 V ze 
活塞 式 驱 痊 按 天 


EA 


(13) 
也 是 CPVCM 的 反例 。 在 孔 孙 尺 度 上 ， 多 相 渗 流 按 活 塞 式 驱 普 的 方式 流动 。 为 了 维护 
CPVCMI5. 1H， 一 些 研究 者 提 0 人 饱和 度 流体 的 属于 未 定义 流体 ， 所 以 它们 的 压力 和 渗流 速度 不 存在 ， 所 以 活 
E 葵 不 能 否定 CPVCM。 这 种 解释 实际 上 忽视 了 常规 的 连续 介质 渗流 理论 采用 了 0 饱和 度 、0 流量 的 事 
实 ; 
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2.2 质疑 2: 在 均 质 储 层 中 CPVCM 条 件 与 Rankine-Hugoniot 条 件 的 矛盾 


两 相 驱 蔡 的 渗流 问题 ， 属 于 流体 力学 的 一 类 激 波 问题 。Rankine-Hugoniot( 跳 跃 ) 条 件 (公式 (14)) Pele 
究 流 体力 学 激 波 问题 的 经 典 理论 公式 ， 不 妨 简 称 R-H 条 件 。 


Uy, — Uy, 
U, 三 Sas (14) 


AF: 5, 表示 含水 体积 比 U, 表示 前 缘 界 面 的 移动 速度 ; u, 表示 流体 (真实 流速 ， 与 渗流 速 有 关系 


WU, 二 Vi /人 办。 根据 公式 (8-9) 得 到 公式 (15) 四 


1 
= glAn -0.5K (8.7. +408) | (15) 
而 按照 CPVCM 给 定 的 界面 条 件 有 : [u -u | N=0 (10) 
ey u, +u a) 
u = 
3 


显然 表达 式 (17) 与 (14) 不 一 致 ， 说 明 CPVCM 与 R-H 条 件 不 相符 合 。 
与 CPVCM 不 同 ，R-H 条 件 可 以 描述 毛细 管 或 孔隙 尺度 下 的 两 相 活塞 式 驱 老 
模拟 中 广泛 应 用 06 17, 19, 24, 29-351, 


2.3 质疑 3: 两 块 毛 管 力 不 同 的 岩 块 中 CPVCM 条 件 与 流体 饱和 度 场 的 矛盾 

岩石 物性 变化 引起 的 突变 界面 属于 固定 界面 ， 但 因为 毛管 力 曲线 不 同 ， 在 CPVCM 条 件 限制 下 ， 

p (S55)= pi (Si )， 所 以 流体 饱和 度 必 须 间断 分 布 。Duijn、Molenaar 和 Neef 建 立 了 两 块 兰 石 毛管 力 不 同 

的 突 主 界面 条 件 & 39， 以 图 5 为 例 ， 有 如 下 结论 : CL) 突变 界面 处 各 相 流体 的 流量 续 ， 总 流量 连续 ，(2) 

湿 相 压力 多 连续 ，(3) 在 在 两 块 岩 块 界面 处 ， 岩 块 的 毛管 力 值 满足 如 下 关系 : 
Aa if (5S- >5") 


S+=0 if (5° <S") 


$ 
w 
一 出 


而 的 移动 ， 在 微观 渗流 


(18) 


Van Duijn 等 人 还 讨论 了 多 种 类 型 的 毛管 力 曲 线 组 合 〈 如 图 6)， 在 此 我 们 将 这 组 突变 条 件 命 名 为 Duijn- 
Molenaar-Neef 突变 界面 条 件 ， 简 称 DMN 界面 条 件 。 这 些 组 合 的 区 别 点 在 上 述 3) AA [AID 3638], DMN 
界面 条 件 在 流体 平衡 状态 下 可 以 通过 实验 室 得 到 验证 ， 也 适用 更 多 的 毛管 曲线 组 合 类 型 (如 图 7)， 但 是 包括 
我 们 的 一 些 研究 者 发 现 它们 不 适用 于 描述 流体 在 非 平稳 状态 下 的 分 布 21, 汶 。 当 通过 技术 手段 控制 流体 出 现 
短暂 的 非 平衡 分 布 。 方式 (1): 首先 让 两 种 岩心 隔离 放置 ， 任 意 配 置 每 块 岩石 让 流体 饱和 度 
Sp Sh €[Syi1-S,,] Sa RAK: Sp 残余 油 )， 但 确保 pi (S7) po (Si); 然后 足够 迅速 地 将 两 块 岩 
心 接触 ， 就 能 突破 DMN 界面 条 件 对 流体 分 布 模 式 的 限定 ， 方式 (2) 利 用 一 些 流体 的 表面 张力 受 磁场 影响 的 
性 质 ， 利 用 磁场 瞬间 完成 毛管 力 曲线 的 改变 ， 从 而 让 流体 由 平衡 状态 瞬间 转 入 非 平衡 状态 。 
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图 5 两 种 不 同 岩 石 的 5 种 毛管 力 曲 线 组 合 类 型 
Fig. 6 Four types of capillary pressure curve compose of two different rocks studied by Van Dujn et al 

在 Van Duijn ÆA (1995 年 ) 的 工作 之 后 ， 所 建立 的 界面 条 件 就 被 用 于 油 藏 数值 模拟 技术 ， 包 括 有 限 
体积 法 和 混合 有 限 元 法 ， 以 提高 毛细 管 压力 曲线 不 同 条 件 下 两 种 岩石 界面 处 的 流量 精度 。 研 究 发 现 通 过 3 
入 DWN 界面 条 件 ， 执 行 毛细 管 平衡 ， 用 一 个 局 部 非 线 性 系统 实现 了 跨 界面 流量 的 计算 90-4， 提高 了 油 藏 
数值 模拟 在 精度 6 和 非 线 性 收敛 方面 的 表现 9 天。 表面 上 这 些 进步 当归 功 于 DMNI, AAPA LAHE 
DWN 和 CPVCM 的 合理 性 ， 但 真相 并 非 如 此 ， 我 们 调查 文献 发 现 引 入 总 速度 -分 数 流 公式 (v, -fy FEA 
续 界面 问题 转换 为 连续 界面 问题 才 是 真正 原因 B91。 

2.4 质疑 4: CPVCM 条 件 不 能 用 于 跨 突变 再 面 的 流量 计算 

油 藏 工程 有 很 多 问题 需要 计算 跨 突变 界面 的 流量 ， 因 为 界面 突变 ， 需 要 采用 数值 方法 近似 计算 。 比 如 
油 藏 数值 模拟 构建 网 格 界面 处 的 数值 流量 、 动 态 法 测试 相 渗 曲线 、 示 踊 剂 、 多 重 介质 和 孔 隐 定 流量 计算 等 等 。 


kais Kish Pas A 
Si» Poi 


G Xan Xu x 
图 6 突变 界面 两 侧 的 微 元 及 它们 的 主要 参数 
Fig. 9 Infinitesimal elements on both sides of the mutation interface and their main parameters 


以 图 9 为 例 ， 并 根据 CPVCM 界面 条 件 〈 见 公式 (4)) 跨 突 变 界面 处 的 水 流量 为 : 


下 Pwi — Prin 
MRTE Ct a ra) (19) 
Ai Ki Aviat AKA; 

当 微 元 计 1 中 水 相 饱 和 度 5,,, (RF SK I AES) ARE =0， 将 其 代入 公式 (19) 有 
gq, 0 ， 说 明 水 无 法 越过 位 置 区 ， 这 与 驱 替 过 程 的 物理 规律 不 相符 合 。CPVCM 推导 出 的 数值 流量 格式 本 
质 上 是 传导 率 系 数 的 调和 平均 ， 在 油 藏 数值 拟 中 被 证 明 是 一 种 失效 的 方法 5， 所 以 被 单 点 上 游 权 (SPU) 
取代 ，SPU 是 目前 商业 油 藏 数值 模拟 软件 采用 的 主要 方法 。SPU 虽然 比 CPVCM 更 合理 ， 但 在 模拟 突变 界 
面 渗流 问题 时 ， 产 生 很 强 的 数值 弥散 ， 误 差 能 够 超过 30% 以 上 5563， 为 此 人 们 又 研究 了 高 阶 精度 方法 ， 包 
ta EDF IODA, Godunovi, TVD, ENO, WENO!®?.81 等 ， 虽 有 改进 但 改进 有 限 ， 并 非 计算 突变 界 
面 流 量 的 理想 方法 。 因 为 它们 并 非 针 对 间断 界面 条 件 的 严格 推导 ， 而 是 不 得 已 的 选择 691。 

用 SPU 和 高 精度 计算 格式 结合 了 总 速度 Vv 公式 -分 数 流 人 公式 避 开 了 间断 界面 问题 50， 成 为 当前 油 
藏 数值 模拟 技术 发 展 的 主流 方法 。 这 促使 我 们 思考 用 总 流速 w 连续 代替 CPVCM 作为 突变 界面 的 处 的 连接 
Kiko HRH v,- 有 ,公式 改造 从 理论 上 只 能 用 于 仅 限 于 不 可 压缩 流体 ， 因 为 对 于 不 可 压缩 流体 理论 上 尚未 
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确定 在 突变 界面 处 的 连续 性 。 
3. 传 统 突 变 界 面条 件 的 推理 的 回 漳 

质量 守恒 定律 被 认为 是 得 出 CPVCM 的 基础 ， 调 研 文献 发 现 渗流 理论 文献 由 质量 守恒 定律 得 到 CPVCM 
有 三 种 推理 方式 : 

(1) 方式 1: 根据 质量 守恒 定律 ， 在 流动 空间 取 任 何 包括 有 突变 界面 的 微 元 〈 见 图 1 所 示 )， 有 如 下 
的 质量 守恒 方程 : 


Min Mirout = Miche +m (20) 


(src 


AP Mg 、1 分 别 表示 表示 单位 时 间 内 流入 和 流出 微 元 中 流体 (的 流量 ，7meng 表 示 单 位 时 间 内 微 元 中 的 
流体 {质量 变化 量 ; 1， 表示 单位 时 间 内 微 元 中 由 源 汇 项 产生 的 流体 ! 质 量 。 传 统 的 渗流 力学 认为 : 当 微 
元 为 面 ， 并 且 不 考虑 源 汇 项 ， 则 公式 (20) 右 侧 得 0， 于 是 有 : 
Mein = Mig OF Pin = PL tout (21) 
对 于 不 可 压缩 流体 ， 密 度 不 变 ， 所 以 有 : 
liin = Vou OF (AV in 1) + (ALY iou Ny) =0 (22) 
所 以 有 [v,]=0 or v, =v; (23) 


但 我 们 认为 ， 上 述 推理 过 程 中 在 利用 质量 守恒 定律 推导 界面 条 件 的 过 程 中 ， 发 生 了 对 象 的 指 代 错 位 ， 
如 图 10 所 示 : 界面 条 件 需要 确定 界面 两 侧 不 同 的 两 个 质点 M _ 、M ”的 速度 关系 ， 而 不 是 突变 界面 上 相 
同 的 质点 M 的 自己 与 自己 比较 速度 ， 即 不 能 将 ” 一石，X Ar yy =A =O 或 者 M =M=M*, 
不 能 由 下 在 一 并 后 三 0 得 由 YW 一 Wr =O ( 即 [”]=0 )， 在 推理 过 程 中 需要 始终 保守 所 比较 的 是 两 个 
不 同 质点 的 参数 ， 对 于 不 能 将 有 厚度 的 三 维 微 元 极限 推理 到 二 维 ，( 或 对 有 厚度 的 二 维 微 元 用 极限 降 维 到 
一 维 )， 否 则 就 是 预 设 了 ” WYP 二 ， 再 证 明 “在 界面 处 连续 就 是 自 证 。 


一 十 
XN 


一 F 
X Xp x 


图 7. 突 变 界 面 两 侧 的 两 个 不 同 质点 与 界面 上 的 质点 的 区 别 B9, 的 ]; 


Fig.10 The two particles on both sides of the sharp interface and the one at the sharp interface’? 64 
(2) 方式 2: 当 流 体 通 过 没有 源 汇 项 的 突变 界面 从 一 个 流动 区 域 Q2 流向 另 一 个 区 域 42 COL 1 所 
示 )， 根 据 物 质 不 灭 原理 ， 从 5 流出 的 量 总 是 等 于 流入 区 域 42, 中 的 流量 ， 所 以 得 出 [ow] = 0 或 [v | =0。 
我 们 认为 方式 “2)〉 的 推理 过 程 同 样 发生 了 对 象 的 指 代 错位 ， 用 Q 区 域 在 界面 工 上 的 质点 M 的 流速 代 
SYM 的 流速 ， 同 样 用 Q, 边界 上 工 的 质点 M 的 渗流 速度 代 蔡 了 M 的 流速 。 根 据 质量 守恒 定律 ， 流 体 
经 过 Q MO, 的 共同 边界 工 从 Q 流出 的 流量 iio 确实 必须 等 于 流入 Q, 区 域 的 流量 TA ， 但 这 规律 只 要 
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D 上 任意 的 流速 ww 可 实现 ， 并 不 要 求 y- =v =v ， 如 果 基于 前 者 就 得 出 后 者 ， 实 质 上 同样 预 置 了 vy, 在 
突变 界面 工 处 连续 ， 然 后 再 自 证 w, 连续 ， 也 是 自 证 。 


> 
o 12 3 4 5 6 7 8 9 10 
时 间 min 


图 8. 突 变 界面 可 视 作 每 个 子 区 域 的 外 边界 的 示意 图 


Fig.8 The inner sharp interface can be regarded as the outer boundary of each sub region 
进一步 以 图 8 为 例 说 明 : 突变 界面 工 既 是 @ 的 外 边界 也 是 区 域 Q, 的 外 边界 。 假 设 在 t<Smin 时 边界 按 
quo, =60cm3/min, gq’, =40 cm3/min 供 液 让 Q, 区 域 保持 与 边界 相同 的 稳 态 流动 ， 然 后 在 t=5min 时 ， 改 变 产 
ÆN q'o =20em3/min, gi, =80 cm3/min， 注 入 条 件 的 改变 将 在 Q, 产生 驱 替 前 缘 ， 前 缘 位 置 为 时 间 的 函数 ， 
用 xy (1) 表示， 可 以 证 明 对 于 任意 的 a ->0 且 e > 0， 总 存在 足够 短 时 间 5， 在 1=5+O(min) 时 刻 ， 满 
Ex, (+6)< 各 +E =x*， 即 能 够 确保 Avi: 而 另 一 端 Q, 区 域 作 为 下 的 上 游 ， 按 工 供 液 的 来 源 ， 只 
需 令 其 按 qr， 和 gi, 相同 的 流量 供 液 ， 根 据 上 游 权 这 一 物理 规律 ， 有 Vi =v HAv Av. 


(3) 方式 3: 以 Hassanizadeh 和 Gray 为 代表 的 理论 证 明 方 式 

以 Hassanizadeh 和 Gray 等 人 的 为 代表 的 学 者 从 理论 的 角度 对 CPVCM 进行 了 证 明 ， 并 受到 广泛 引用 2 
4,65]。 与 方式 《1-2) 的 推理 方式 不 同 ， 他 们 的 推导 过 程 中 严格 遵守 了 x 和 xx 指向 不 同 的 质点 。 但 我 们 发 
现在 其 证 明 过 程 中 没有 考虑 到 微 元 厚度 b 改变 《垂直 于 突变 界面 方向 ， 如 图 1 所 示 ) 将 影响 微 元 的 平均 饱 
AIRES,» BUDS, 应 为 非 0 常数 ， 所 以 不 能 得 出 : 


mm 


_ a 
lim (b (P), 95, )=0 (24) 


进而 不 能 得 出 V =V o PRANAB SR SCARE 6112, DAZ Peng EA (2019) 对 其 做 的 核对 B9 。 
总 之 通过 文献 追 足 和 分 析 ， 我 们 认为 传统 渗流 力学 在 用 质量 守恒 方程 推导 突变 界面 条 件 中 将 界面 上 的 
党 速 假设 为 界面 两 侧 的 流速 ， 通 过 自 证 的 方式 建立 了 CPVCM， 是 不 合理 的 。 
4. 突变 界面 条 件 的 重新 建立 
本 节 我 们 仍然 根据 质量 守恒 定律 和 两 相 流 达 西 公式 建立 突变 界面 条 件 。 在 前 期 我 们 针对 一 维 问题 开展 
THRO 的 ， 在 此 进一步 扩展 到 三 维 。 
4.1 跨 突 变 开 面 的 总 流量 及 各 相 流量 计算 
如 图 1 所 示 ， 油 藏 被 突变 界面 工分 割 为 两 个 连续 子 区 域 D 、Q,。 界 面 P 、T 在 工 两 侧 并 与 工 平行 
而 且 无 限 靠近 I 。 在 突变 界面 工 两 侧 取 两 点 X“ (x,y,z) MXA (x,y,z) ， 两 点 连 线 记 为 1.p ， 


< 


X (x,y,z) A lg 与 突变 界面 工 的 交点 ，X (x,y) Alg ST MBE, Xt (x,y,z) Blog 
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ST HZA, WREX >X", XX”. ERX MX? HERE: 


BEET, 长度 大 于 线段 -+ ， 


毛管 力 和 压力 梯度 足够 相近 
(5) 得 到 两 相 总 流速 运动 方 和 


v, =v, +v, =-KA (Vp, -上 Vp -万 gVD) 
引入 全 拟 压力 (global pressure) 的 概念 [66. 67 


1 6p.(5) 
=p + 一 一 一 

Pn f fy (5) ôç 
将 〈26) 代入 等 式 25) 得 并 对 等 式 两 XP RIRS: 


x’ 
| 


dç 


| )-al 

“SKA, 
根据 ,无 散场 的 连续 性 质 ，X“ 入 4 ERRET, lag 线段 上 V 近似 相等 ， 
置 于 积分 符号 外 ， 用 nag 表示 沿 1. 的 单 人 


位 矢量 。 所 以 得 到 : 
D -9° + pg(D’ -D*) 
(mojF r 
i —~K°)° +]. —K* AP 
HF Lg 和 界面 工 的 法 线 方向 一 致 ， 所 以 ng 三 让 所 以 跨 突变 界面 的 总 流 


量 


all 


qr =v =(v, nr )Ar = A | ©” -@* +g (D° -D*) | 


x dl [rw dl 
lee | 


K°Ae +x" —KPAP 
HAR 28) 可 计算 在 X M Xt 处 各 相 渗流 沿 工 法 


线 方向 的 渗流 速度 : 


K = f2v,-GP* v= ftv, +GP" 


vi = Jp -GP? vt = fv, +GP’ 


比较 公式 30) 中 4 项 ， 可 发 现 一 


般 情 况 下 突变 界面 两 侧 各 相 渗流 速度 并 不 相等 ; 
v EV, V +V 


如 何 确 定 突变 界面 工 上 各 相 渗流 速度 则 根据 流体 力学 的 上 游 权 法 则 


， 所 以 有 : 
i =upstream(v,.v;) =if (® > ®*;vi;v} ) Ce{n,w} 
: b=0+pgaD=[ 人 (5 jle) ) 


dc+ + AD 
| ad (fp,+f.p,)8 
所 以 跨 突 变 界 面 的 各 相 流 量 为 : 


其 


qn, =y, “nA, ae Daa —Aro, „OPZ; fr qr 


-Ar0,,GP’) 
q =v; nA = 六 (> 8; fq 


T +A.0,, GP: fig! +4,0,, GPP) 
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X” RX? AGES XT 相隔 太 远 ， 以 确保 线段 EX” AEX Y 


BBL, 与 工法 


线 方向 相同 ; 


mi. AE 


根据 (25-a) 得 到 Vv, =0, Bl vy, 在 整个 渗流 场 为 无 散场 或 保守 场 。 由 公式 


(25) 


(26) 


(27) 


所 以 等 式 (27) 左 侧 y 可 


(28) 


(29) 


(29) 


(30) 


(31) 


(32) 


(33) 
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3 跨 突变 邢 面 的 处 压力 的 分 布 


根据 等 式 (28) 很 容易 证 明 全 拟 压 力 (global pressure) 在 突 界面 处 连续 。 即 : 
lim © = 中 ! = lim ®* (34) 


说 明 由 等 式 (26) 所 定义 的 全 拟 压力 (global pressure) 在 突 界面 处 连续 。 
DAHA XM XHA XP X" 两 点 按 (28) 式 计算 总 流速 在 线段 1 方向 分 量 (v, Nap) 
并 联 立 求解 可 工 处 的 全 拟 压力 中 : 


x? dl 
p fo K’ Af 


o o (®* + pg (D° -D")) 
上 K°Ar he Kia 


(35) 


"dl 
CS 


ee ee 
x dl x” dl (0° +pe(D -D")) 
i K° A? a KAP 
根据 〈36) MWR AWE MATOS, HTTP SEL, 线段 上 和 AX Wb HE) 


ove ape (5) ap: (5) 
ôç 


=® — r sele) de m=" +f” f(e) dç 
(36) 


‘ . 1 4, êp 
本 £2 (6) a p=" +] EA 


根据 (36) 式 可 知 ， 当 突 界面 两 侧 f,、 p 等 参数 不 同 ， 得 到 的 p,, 和 p, 在 界 突 主 面 处 不 连续 。 结合 
和 (36) 形成 的 突变 界面 处 各 相 压 力 和 渗流 速 不 连续 ， 由 此 构成 的 界面 条 件 不 同 于 CPVCM， 为 示 区 
其 称 为 卫 VCM。 
突变 界面 及 邻近 区 域 渗流 速度 的 计算 


公式 (29) 
别 ， 将 


公式 (28) We ,方向 的 总 流速 分 量 。 对 于 一 维 渗流 问题 ， 该 分 量 即 为 总 流速 分 量 ， 
但 对 于 二 维和 三 维 渗流 问题 ， 还 需要 确定 其 它 方向 的 总 流速 分 量 。 Tmax (x y”, z” ) 作 为 局 部 坐标 


原点 ， 在 界 ee yeh), XX GT RSL, BOMEX" (x",y", 


z’) 


XE (2%, y%, 2°), HAERERE: XIXE GOAL), XIX” GOW.) X°X Aly) 在 空间 
相互 垂直 ， 如 果 他 们 与 工 有 交点 ， 记 为 X2(B s{ 生 75) WRH, X’ (eléng) 5X? 位 于 


同一 的 子 区 域 。 


中 文 版 阶段 的 稿件 


图 9 在 三 维 渗 流 空间 计 香 


总 流速 张 量 的 示意 


图 


Fig. 13 schematic diagram of calculating total velocity in three-dimensional porous media 


用 三 次 公式 〈38) 计算 三 个 方向 的 总 流速 的 分 量 : 


D“ -9° + pg( D" -D*) 


(v, Nag )= xT 1 


其 中 e(X5)=if (3X710). 


xh 


一 < 
(Xe) fye gage Ka) fs gage (tel Kea) fe KAP 


—BAEMSY on 


再 由 三 个 分 量 组 成 总 流速 的 张 量 矩阵 (39) : 


ve £ v > ve vn” Vern ven 
Bay lal Ve | ae ane (38) 
v v% i ve v&n” ya v n” 
B_ ya Pp_ ,a Bia 
FH pe =~ * , ye YY, yh at (Jelent) 
pea] ex]! fo 
进一步 根据 公式 (89) 可 以 计算 和 处 各 相 流速 。 将 (38) (38) FRAZS zh(29) Ay EARMA XTA 


相 流速 。 
4.5 渗流 方程 中 突变 宽 面 条 件 的 重新 建立 


AST HES JPVCM 采用 了 : (1) 总 流速 连续 VY.v, =O; (2) 界面 处 的 流速 取 上 游 权 。 综 合 考虑 CPVCM TF 
在 的 问题 和 JPVCM 的 建立 过 程 ， 我 们 在 多 相 渗流 方程 中 ， 连 接 两 个 连续 子 区 域 的 突 主 界面 采用 如 下 的 界 
面条 件 〈 见 公式 40): 

vl- =if (® >"; v; vi) (a) Bo 


v +v, =v +v (b) 
理由 如 下 : 
C1) 所 列 界 


(2) 这 是 两 项 被 渗流 力学 认可 的 结论 : Vv =0 根 据 质量 守恒 得 到 ; 


用 条 件 可 以 确定 界面 两 侧 流 体 压力 和 渗流 速度 等 参数 ， 实 现 对 相 邻 子 区 域 的 渗流 方程 的 连接 ; 


EF 游 权 在 流体 力学 中 是 一 种 自然 规 


律 的 表达 ， 在 建立 渗流 数学 模型 前 ， 它 就 应 该 列 入 物理 模型 ， 
有 对 应 的 数学 公式 表达 ， 而 非 作 为 一 种 可 选 的 数值 求解 方法 ; 


(3) 公式 (40-a) 根据 中 确定 上 游 权 方向 。 相 比 一 些 文献 采用 


作为 前 提 假 设 ， 在 所 建立 的 渗流 数学 模型 中 


P: = Pi 十 P18AD 确定 上 游 权 ， 它 更 为 合 
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理 ， 因 为 当 两 相 流 体 对 流 时 ， 利 用 p, (40-a) 判断 上 游 得 出 的 结果 可 能 违背 界面 条 件 (40-b); 采用 有 限 
小 时 间 间 隔 分 析 法 ， 在 间隔 内 ~ lim (f +Af) 时 间 内 ， 用 中 判定 的 上 游 函 数 确定 的 y, ，v, 变化 平缓 ， 


但 是 由 五 ,判定 上 游 权 结果 会 出 现 跳跃 ， 如 图 2、 图 4、 图 9、 图 11 所 示 。 
S. 结 论 


(1) 传统 渗流 力学 理论 认为 根据 质量 守恒 定律 在 突变 界面 处 各 相 渗流 速度 连续 (CPVCM)， 本 文 从 三 
个 方向 进行 了 质疑 : 根据 CPVCM 推断 的 突变 界面 两 侧 的 流体 饱和 度 分 布 存在 实际 的 反例 ， 包 括 均 质 介质 
和 毛管 力 不 同 的 两 种 多 孔 介质 两 类 ; 对 于 驱 蔡 前 缘 问 题 CPVCM 与 Rankine-Hugoniot 跃 变 条 件 矛 盾 ; 
CPVCM 推导 出 的 数值 计算 格式 被 油 藏 数值 模拟 证 明 会 出 现 非 物理 结果 而 被 弃 用 ; 

(2) 论文 回溯 了 文献 根据 质量 守恒 方程 推导 CPVCM 的 过 程 ， 发 现 这 些 推 导 过 程 出 现 用 界面 上 同一 质 
点 的 流速 代替 了 界面 两 侧 不 同 质点 的 流速 ， 实 质 上 预 置 了 渗流 速度 连续 ， 属 于 自 证 。 
(3) 本 节 根 据 质量 守恒 方程 、 达 西 公式 、 上 游 权 性 质 ， 针 对 不 可 压缩 流体 ， 考 虑 两 相 流 毛 管 力 和 重 
力作 用 ， 在 三 维 空间 建立 了 突变 界面 两 侧 各 相 流速 ， 得 出 在 突变 界面 两 侧 ， 各 相 渗流 速度 和 各 相 压 力 一 般 
都 不 连续 ， 称 为 了 了 VCM。JPVCM 能 够 解释 CPVCME 面 临 的 上 述 所 有 质疑 ， 不 仅 实 用 于 有 间断 界面 的 渗流 ， 
也 适用 于 连续 场 渗流 ， 而 且 也 适用 于 纯 活 塞 式 驱 蔡 ; 

(4) 根据 JPVCM 突变 界面 处 的 变 界 面条 件 包 括 两 项 : 其 一 、 突 变 界面 处 总 流速 连续 ;其 二 、 突 变 界 
面 处 的 渗流 速度 由 上 游 权 函 数 确 定 。 
当前 关于 多 相 渗流 的 突变 界面 的 条 件 的 研究 仅 限 于 不 可 压缩 两 相 流 ， 而 考虑 流体 可 压缩 性 、 相 间 传 质 
或 相 变 、 三 相 流 、 黑 油 模型 、 组 分 油 藏 模型 等 更 复杂 情况 还 有 待 研 究 。 

实际 油 藏 基本 含有 变 变 界面 ， 油 藏 工程 方法 必须 依赖 考虑 突变 界面 的 渗流 力学 ， 所 以 突变 界面 条 件 完 
善 作为 油 藏 工程 的 基础 理论 ， 将 影响 油 藏 开发 的 各 个 环节 ， 例 如 相对 渗透 率 测量 、 试 油 试 井 、 油 藏 数值 模 
拟 、 生 产 动态 分 析 、 产 量 预 测 、 开 发 方案 优化 、 油 气 二 次 运 移 成 藏 、 井 简 计 算 、 微 观 渗流 等 。 除 了 油气 层 
渗流 ， 也 将 在 其 它 多 相 渗流 领域 发 挥 作 用 。 


T 


s 


Hy 


X 


参考 文献 
[1] MUSKAT M, WICKOFF R D. The Flow of Heterogeneous Fluids Through Porous Media [J]. The Journal of Geology, 
1937, 46(2): 


[2] SHELDON J W, ZONDEK B, CARDWELL W T. One-Dimensional, Incompressible, Noncapillary, Two-Phase Fluid 
Flow in a Porous Medium [J]. Transactions of the AIME, 1959, 207(1): 136-43. 


[3] BACHMATY B, J. Mathematical formulation of transport phenomena in porous media; proceedings of the the Second 
Symposium on Fundamentals of Transport Phenomena in Porous Media, IAHR, Guelph, Ont., Canada, F, 1972 [C]. 


[4] BEAR J. Hydraulics of Groundwater [M]. New York: McGraw-Hill, 1979. 


[5] HASSANIZADEH M, GRAY W G. General conservation equations for multi-phase systems: 1. Averaging procedure 
[J]. Advances in Water Resources, 1979, 3(1): 25-40. 


[6] HASSANIZADEH S M, GRAY W G. Derivation of conditions describing transport across zones of reduced dynamics 
within multiphase systems [J]. Water Resources Research, 1989, 25(3): 529-39. 

[7] RAATS P. The role of inertia in the hydrodynamics of porous media [J]. Archive for Rational Mechanics & Analysis, 
1972, 44(4): 267-80. 

[8] BEAR J. Dynamics of Fluids in Porous Media [M]. Courier Corporation, 1975. 


[9] VAN DUIJN C J, MOLENAAR J, DE NEEF M J. The effect of capillary forces on immiscible two-phase flow in 
heterogeneous porous media [J]. Transport in porous media, 1995, 21(1): 71-93. 


12 


中 文 版 阶段 的 稿件 


[10] 刘 绍 学 , 朱 元 森 . 数学 辞海 [M]. 数学 辞海 , 2002. 
[11] HASSANIZADEH S M, GRAY W G. Boundary and interface conditions in porous media [J]. Water Resources 
Research, 1989, 25(7): 1705-15. 


[12] HASSANIZADEH S M, GRAY W G. Mechanics and thermodynamics of multiphase flow in porous media including 
interphase boundaries [J]. Advances in Water Resources, 1990, 13(4): 169-86. 


[13] JAFFRé J. Numerical calculation of the flux across an interface between two rock types of a porous medium for a 
two-phase flow [M]. Hyperbolic Problems Theory Numerics Applications. NY; Stony Brook. 1996: 165-77. 


[14] LIU Z. Propagation and Evolution of Wave Fronts in Two-Phase Porous Media [J]. 1999, 34(1): 209-25. 

[15] FUCIK R, MIKYSKA J, BENES M, et al. Semianalytical Solution for Two-Phase Flow in Porous Media with a 
Discontinuity [J]. Vadose Zone Journal, 2008, 7(3): 1001-7. 

[16] ALI, Q., RAEINI, et al. Modelling two-phase flow in porous media at the pore scale using the volume-of-fluid method 
[J]. Journal of Computational Physics, 2012, 231(17): 5653-68. 

[17] RAEINI A Q, BIJELJIC B, BLUNT M J. Generalized network modelling: capillary-dominated two-phase flow - 
model description [J]. PHYSICAL REVIEW E, 2017, 97(2): 023308. 

[18] DEDè L, QUARTERONIA. Isogeometric Analysis of a Phase Field Model for Darcy Flows with Discontinuous Data 
[J]. 数学 年 刊 : B 辑 英文 版 , 2018, 39(3): 26. 

[19] SHAMS M, RAEINI A Q, BLUNT M J, et al. A numerical model of two-phase flow at the micro-scale using the 
volume-of-fluid method [J]. Journal of Computational Physics, 2018, 357(159-82. 

[20] ALALI AH, HAMON F P, MALLISON B T, et al. Finite-volume simulation of capillary-dominated flow in matrix- 
fracture systems using interface conditions [J]. Computational Geosciences, 2021, 25(1): 17-33. 

[21] TRAN LK, KIM J C, MATTHIS K. Simulation of two-phase flow in porous media with sharp material discontinuities 
[J]. Advances in Water Resources, 2020, 142(1): 103636. 

[22] DASHTBESH N, ENCHéRY G, NOETINGER B. A dynamic coarsening approach to immiscible multiphase flows 
in heterogeneous porous media [J]. Journal of Petroleum Science and Engineering, 2021, 201(3): 108396. 

[23] BOSMAS B, HAMON F P, MALLISON BT, et al. Smooth implicit hybrid upwinding for compositional multiphase 
flow in porous media [J]. Computer Methods in Applied Mechanics and Engineering, 2022, 388(114288. 

[24] SHOKRI J, GODINEZ - BRIZUELA O E, ERFANI H, et al. Impact of displacement direction relative to 
heterogeneity on averaged capillary pressure - saturation curves [J]. Water Resources Research, 2022, e2021WR030748. 


[25] XIN J, SHI F, FAN S, et al. A sharp interface multiphase flow model for two-dimensional water impact of a symmetric 
and asymmetric wedge [J]. Applied Ocean Research, 2022, 119(102988-. 

[26] BUCKLEY S E, LEVERETT M C. Mechanism of Fluid Displacement in Sands [J]. Transactions of the AIME, 1942, 
146(01): 107-16. 

[27] DOSTER F, HILFER R. Generalized Buckley—Leverett theory for two-phase flow in porous media [J]. New Journal 
of Physics, 2011, 13(12): 123030. 

[28] RANKINE W. On the thermodynamic theory of waves of finite longitudinal disturbance [J]. Philosophical 
Transactions of the Royal Society of London, 1870, 160): 277-88. 

[29] JAFARI I, ROKHFOROUZ M R. Numerical modeling of water oil two-phase flow during counter-current 
spontaneous imbibition in porous media at pore-scale [J]. Petroleum Science and Technology, 2020, 38(24): 1040-53. 
[30] WANG D, CHENG L, CAO R, et al. Pore-scale model of two phase flow in 2D porous media: Influences of interfacial 
tension and heterogeneity effects on CO2 injection in the tight oil reservoir [J]. IOP Conference Series: Earth and 
Environmental Science, 2020, 467(1): 012061. 

[31] NOIRIEL C, SOULAINE C. Pore-Scale Imaging and Modelling of Reactive Flow in Evolving Porous Media: 
Tracking the Dynamics of the Fluid—Rock Interface [J]. Transport in Porous Media, 2021, 140(1): 181-213. 

[32] RAEINI Q, ALI. Modelling multiphase flow through micro-CT images of the pore space [D]; Imperial College 
London, 2013. 

[33] RAEINI A Q, BLUNT M J, BIJELJIC B. Direct simulations of two-phase flow on micro-CT images of porous media 
and upscaling of pore-scale forces [J]. Advances in Water Resources, 2014, 74(116-26. 

[34] MEHMANI Y, TCHELEPI H A. Multiscale Formulation of Two-Phase Flow at the Pore Scale [J]. Journal of 
Computational Physics, 2019, 389(164-88. 

[35] SAHAR, BAKHSHIAN, SEYYED, et al. Pore-scale characteristics of multiphase flow in heterogeneous porous 
media using the lattice Boltzmann method [J]. Scientific Reports, 2019, 9(1): 1-13. 

[36] NEEF M, MOLENAAR J. Analysis of DNAPL infiltration in a medium with a low-permeable lens [J]. Computational 
Geosciences, 1997, 1(2): 191-214. 


中 文 版 阶段 的 稿件 


[37] DUIJN C, MOLENAAR J, NEEF M. The effect of capillary forces on immiscible two-phase flow in heterogeneous 
porous media [J]. Transport in Porous Media, 1995, 21(1): 71-93. 

[38] VAN DUIJN C J, DE NEEF M J. The Effect of Capillary Forces on Immiscible Two-phase Flow in Strongly 
Heterogeneous Porous Media [J]. 1994, 

[39] PENG X, MO F, LIANG B, et al. Critical Investigation Of Interface Conditions For Fluid Pressures, Capillary 
Pressure, And Velocities At Jump Interface In Porous Media [J]. Journal of Porous Media, 2019, 22(6): 723-44. 

[40] AHMED R X Y, EDWARDS M G. A cell-centred CVD-MPFA finite volume method for two-phase fluid flow 
problems with capillary heterogeneity and discontinuity [J]. Transport in Porous Media, 2019, 127(1): 35-52. 

[41] BRENNER K, CANCéS C, HILHORST D. Finite volume approximation for an immiscible two-phase flow in porous 
media with discontinuous capillary pressure [J]. Computational Geosciences, 2013, 17(3): 573-97. 

[42] CANCES C. Finite volume scheme for two-phase flows in heterogeneous porous media involving capillary pressure 
discontinuities [J]. Esaim Mathematical Modelling & Numerical Analysis, 2009, 43(5): 973-1001. 

[43] ENCHéRY G, MICHEL R E. Numerical Approximation of a Two-phase Flow Problem in a Porous Medium with 
Discontinuous Capillary Forces [J]. Siam Journal on Numerical Analysis, 2006, 43(6): 2402-22. 

[44] FRIIS, H. A, EVJE, et al. Numerical treatment of two-phase flow in capillary heterogeneous porous media by finite- 
volume approximations [J]. International journal of numerical analysis and modeling, 2012, 9(3): 505-28. 

[45] HAMON F P, MALLISON B T, TCHELEPI H A. Implicit Hybrid Upwinding for two-phase flow in heterogeneous 
porous media with buoyancy and capillarity [J]. Computer Methods in Applied Mechanics and Engineering, 2018, 331(701- 
27. 

[46] ALALI A, HAMON F, MALLISON B, et al. Finite-Volume Simulation of Capillary-Dominated Flow in Matrix- 
Fracture Systems using Interface Conditions [J]. 2019, 

[47] BRENNER K, DRONIOU J, ROLAND M, et al. Total-velocity-based finite volume discretization of two-phase Darcy 
flow in highly heterogeneous media with discontinuous capillary pressure [J]. IMA Journal of Numerical Analysis, 2021, 
[48] BRENNER K, ROLAND M, QUENJEL E H. Vertex Approximate Gradient discretization preserving positivity for 
two-phase Darcy flows in heterogeneous porous media [J]. Journal of Computational Physics, 2020, 409(109357-. 

[49] BRENNER K, CHORFI N, MASSON R. Sequential implicit vertex approximate gradient discretization of 
incompressible two-phase Darcy flows with discontinuous capillary pressure [J]. Computational Geosciences, 2021, 26(1): 
147-69. 
[50] JIANG J, TCHELEPI H A. Nonlinear Acceleration of Sequential Fully Implicit (SFI) Method for Coupled Flow and 
Transport in Porous Media [J]. 2019, 352:246-75. 

[51] WATANABE S, LI Z, BRATVEDT K, et al. A Stable Multi-phase Nonlinear Transport Solver with Hybrid Upwind 
Discretization in Multiscale Reservoir Simulator; proceedings of the Ecmor Xv-european Conference on the Mathematics 
of Oil Recovery, F, 2016 [C]. 

[52] AGHILI J, BRENNER K, HENNICKER J, et al. Two-phase Discrete Fracture Matrix models with linear and 
nonlinear transmission conditions [J]. GEM-International Journal on Geomathematics, 2018, 10(1): 1-35. 


[53] BRENNER K, GROZA M, JEANNIN L, et al. Immiscible two-phase Darcy flow model accounting for vanishing and 
discontinuous capillary pressures: application to the flow in fractured porous media [J]. Computational Geoences, 2016, 
21(3): 1-20. 

[54] BRENNER K, HENNICKER J, MASSON R, et al. hybrid-dimensional modelling of two-phase flow through 
fractured porous media with enhanced matrix fracture transmission conditions fracture model multiphase hybrid- 
dimensional darcy flow model vertex approximate gradient scheme polyhedral meshes [J]. Journal of Computational 
Physics, 2018, 357(100-24. 

[55] AZIZ K, SETTARI A. Petroleum reservoir simulation [J]. Elsevier Applied Science Publishers, 1986, 

[56] BRENIERT Y, JAFFRTT J. upstream differencing for multiphase flow in reservoir simulations [J]. SIAM journal on 
numerical analysis, 1991, 28(3): 685-96. 

[57] TADMOR E. Numerical viscosity and the entropy condition for conservative difference schemes [J]. Mathematics of 
Computation, 1984, 43(168): 369-81. 

[58] WARRICK A W. Numerical approximations of Darcian flow through unsaturated soil [J]. Water Resources Research, 
1991, 27(6): 1215-22. 

[59] RAITHBY GD. A critical evaluation of upstream differencing applied to problems involving fluid flow [J]. Computer 
Methods in Applied Mechanics and Engineering, 1976, 9(1): 75-103. 

[60] TODD M R, APOS, DELL PM, et al. Methods for Increased Accuracy in Numerical Reservoir Simulators [J]. Society 
of Petroleum Engineers Journal, 1972, spe-3516(06): 515-30. 

[61] DICKS E M. Higher Order Godunov Black-Oil Simulations for Compressible Flow [D]; University of Reading, 1993. 


14 


中 文 版 阶段 的 稿件 


[62] BRANTSON E T, JU B, WU D, et al. Stochastic porous media modeling and high-resolution schemes for numerical 
simulation of subsurface immiscible fluid flow transport [J]. Acta Geophysica, 2018, 66(3): 243 -66. 

[63] ABOTYAR T, GEORGOULIS E H, ISKE A. Adaptive ADER Methods Using Kernel-Based Polyharmonic Spline 
WENO Reconstruction [J]. Siam Journal on Scientific Computing, 2010, 32(6): 3251-77. 

[64] PENG X, MO F, DU Z. Interface condition for the Darcy velocity at the water-oil flood front in the porous medium 
[J]. PLoS ONE, 2016, 12(5): e0177187. 

[65] BRENNER K, CHORFI N, MASSON R. Sequential implicit vertex approximate gradient discretization of 
incompressible two-phase Darcy flows with discontinuous capillary pressure [J]. Computational Geosciences, 2022, 26(1): 
147-69. 

[66] ANTONCEV S N. The solvability of boundary value problems for degenerate equations of two-phase filtration [J]. 
Dinamika Splon Sredy, 1972, 28-53,246. 

[67] G CHAVENT J J. Mathematical models and finite elements for reservoir simulation:single phase, multiphase, and 
multicomponent flows through porous media [M]. ELsevier, 1986. 


[68] CHEN Z. Degenerate Two-Phase Incompressible Flow: I. Existence, Uniqueness and Regularity of a Weak Solution 
[J]. Journal of Differential Equations, 2001, 171(2): 203-32. 


